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able of contents 
@ 51 HISTORY ; Detailed Current Edit History 
(5) 5 DECLARATIO S$ 
(4) 108 OTSSPOWGJ - G REAL*8 ** INTEGER*4 
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H 16 
OTSSPOWGJ - G REAL*8 ** INTEGER*4 power routine 16-SEP-1984 91: 3 ound 
1-005 6 (CMTHRTL.SR 610 Hy J.MAR;1 (1) 
G 
0 


6-SEP-1984 


59:3 
1:28:1 
E seve *¢ ae | a INTEGER*4 power routin 
T ; File: OTSPOWGJ.MAR Edit: SBL 1005 


| adannadebnensdesensennbeseebanetieetsenerbedsnunesnesonebnonsenesnndenenenes 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, - MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


*® + 
*® a 
® ® 
‘® ® 
** * 
:* THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED + 
ft ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
> LUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER: 
. JES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY + 
ie PERSON. TITLE’ TO AND OWNERSHIP OF THE SOF TUARE 76 HEREBY * 
:® TRANSFERRED. . 
® ® 
;* THE INFORMATION IN THIS SOFTWARE IS SU3JECT TO CHANGE WITHOUT NOTICE * 
;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
te CORPORATION. . 
*® ® 
*® ® 
** ® 
*® ® 
*® ® 
*?® ® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


COOCooooooooooooooooooooooo 
COOQOoQooooooooQoooooooooooooo 


ARR RRRAAARSALALASLALALASELSELE SESE SELES ERE E SEES RRR EERE E RRR RRR E ERE S EEE R ERS OY 


; : FACILITY: Language support library - user callable 
; : "ABSTRACT: 


G REAL*8 base to INTEGER*4 power. 
Flosting overflow and underflow can occur. 
Undefined exponentation can occur if base is 0 and power 
is 0 or negative. 


; ; VERSION: 1 


: HISTORY: 
: ; AUTHOR: 
Steven B. Lionel, 6-Feb-79: Version 1 


SOoOOOCOCOOSOSOOSOSSOSOOOSOSOOSOSSOSOSOSOSOSOOOOSOSOOOOOOOOOOOOOOOOOO 


SOCoooooooooooooooooooao 
CGOOOQCOOOOCSOSOCSOSOSSSSOSSSSOSSSCOCSOOSCOSOSOOSOOOOOOOOOOOOOOOOOO 


WOON AUS WN = OO DNAUES WIN (OO ONAU EWN 90 OONOAOUE WN OOOVNOUS Wh 
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ar tt HISTORY ; Detailed current Edit History ae ae ts 1 01:38:33 MTHRTL.SRCJOTSPOWGJ .MAR; 1 wt (3) 


5 -SBTTL HISTORY ; Detailed Current Edit History 


dit Wistory for Version 1 of OTS$POWGJ 

- § SBL1002 = Declare externals. BL 17- gr, '979 

1- - Correct some comments. JBS 30-JUL-1979 

1-001 = Adapted from OTS$POWDJ version 1-001. SBL 06-Feb-79 
4 = Add handlers to catch SS$_FLTOVF and SS$ 


C _FLTDIV, and erqnet 
MTHS_FLOOVEMAT or MTHS_FLOUNDMAT instead, depending on the context. 
Also "disable 1V and change BLBS/ASHL at EXPGTR to BBSC/ROTL for 
uniformity with OTSSPOWRJ. JAW 26-Feb-1980. 
1-005 = Use general mode addressing. SBL 30-Nov-1981 


TR Te TE TE Te Tere re Ty 
— 
uy 
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16-SEP-1984 01:59:33 VAX/VMS Macro V04-00 Pp , 
greeks P1isaita PATANTS Gace Os OO ans1 28? Gd, 
.SBTTL DECLARATIONS 


: INCLUDE FILES: 


EXTERNAL SYMBOLS: 


-DSABL GBL 

eEXTRN MTHSK en MTHSK -FLOOVEMAT MTHSK_FLOUNDMAT 

-EXTRN MTHSSSIGN Math error routine 

~-EXTRN SS$_ PLTOVE sss FLTOVF fF; ses _FLTDIV, ¢3$. FLTDIV. F, SS$_CONTINUE 
MACROS: 

SCHFDEF 3; Define condition handler symbols. 

SSFDEF : Define stack frame symbo!s. 

SPSLDEF ; Define program status Longword 


symbols. 


EQUATED SYMBOLS: 


base = 4 


base input formal - by-value 
exp = 12 


exponent intpu formal - by-value 

Note: G_floatin a violates 
calling standar ut ok since this 
routine is a code support routine (OTS$) 


OWN STORAGE: 


PSECT DECLARATIONS: 


-PSECT _OTSSCODE PIC,SHR,LONG, EXE, NOWRT 
3; program section for OTS$ code 
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1805 OTSSPOWG) = G REALSB s* GNTEGERSE grSEei19be 1isesia FATUNTS eaesoreccucs.mars1 29° cf, 
a9 8 -SBTTL OTSS$POWGJ - G REAL*8 ** INTEGER*4 
19 +4 
i 12 ; FUNCTIONAL DESCRIPTION: 
0 Io 2 G_ floating result = G_floating base ** signed longword e ent 
4 13 ; The 6 floating result is elven by: . - Secon 
009 1 ; base exponent result 
60 8 3 any >0 product (base * 2**i) where i is each 
B98 }? ; non-zero bit position in exponent 
a ot ed Jsttnad exponentat 
° = s n n exponen 
$00 198: <b 6088 Wa each oe race lege: 
000 4; 
000 a2 >0 <0 1.0 / product (base * 2**j) 
000 6; where i is each non-zero bit position 
000 3 in texponent! 
000 B 3 = 0 <0 Undefined exponentation 
000 9; <0 <0 1.0 / product (base * 2**i) 
000 0; where i is each non-zero bit position 
38 3 : in texponent: 
00 g ; Floating overflow can occur on either of the two MULG's. If this 
00 34; napoens when the exponent is less than zero, the oncept ten is caught by 
00 a local condition handler named EXC_HNDLR_UNDER, which sets the result 


: 
0 
0 
0 
0 
0 
0 
Q 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
" 
0 
0 
0 
0 


; to 0.0 and either signals MTH$_FLOUNDMAT Tif FU is enabled in the 

3 caller's PSW) or continues at POWGJX. If it happens when the exponent 
3 is greater than zero, the exception is cought by a local condition 

; handler named EXC HNOLR_OVER which sets t 

: operand (-0.0) and signals MTHS$_FLOOVEMAT. 

$ Floating overflow and floating divide by zero can occur on the DIVG. 

3 These exceptions are caught by EXC_HNDL Mg the result to 


Undefined exponentiation occurs if base is 0 and 
exponent is 0 or negative. 


; CALLING SEQUENCE: 
Power.wg.v = OTSSPOWGJ (base.rg.v, exponent.rl.v) 

; INPUT PARAMETERS: 

3 base 
exponent 


; IMPLICIT INPUTS: 
; The setting of FU in the caller's PSwW. 


; OUTPUT PARAMETERS: 
; NONE 


e result to the reserved 


the reserved operand (-0.0) and signals ATHS_ 


- G REAL*8 base 
- INTE 


GER*4 exponent 


ce ce ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ed eed eed 


SOOOoQoooooooooooooooooooosooosoeo 


; IMPLICIT OUTPUTS: 
: NONE 
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OTSSPOWGJ - G REAL* NTEGER oreEpm obs $iiaeita EMTARTE. eaedoreoGucs.mar;1 2% 
0 165 ; 
0 196 : FUNCTION VALUE: 
09 168 3 G_floating base ** signed lLongword exponent 
00 176 ; SIDE EFFECTS: 
00 17 ; Stgnats MTHS FLOOVEMAT if theeting overflow occurs on either of the two 
0000 \t? ; ULG" s when ox onent > 0, or if floating overflow or divide by zero 
000 175 ; Signals te a gee iY Moy. pyecties gecurs - either of the two 
464 178 : ULG's when +f ooagt and caller has FU 
0 77 3 SIGNALS MTHS UNDE (82 = ‘, att EXPONENTAT TOW 5 if 
0000 178; base is 0 and an is O or negative. 
0000 179; 
0000 180 ;-- 
0000 181 
0000 183 
001C $609 iB -ENTRY OTSS$POWGJ, “M<R2, R3, R4> 
OOS 185 ; Disable integer overflow. (Occurs on 
000 186 3 maximum negative exponent.) 
C'AF 9E BoE 13 MOVAB B*EXC_HNDLR_OVER, (FP) ; prnesi ese snsquttans to 
08 SOFD 0006 189 MOVG #1, RO : RO/RT = initial result 
04 AC 5OFD QOOA 190 MOVG base(AP), R2 3 R2/R3 = 
0c AC DO QOOF 191 MOVL exp(AP), R4 3 R4 = tethaore 3 
14 0013 136 BGTR EXPGTR $ prened if exponent > 0 
5D'AF 9E& 0015 19 MOVAB B*EXC_HNDLR_UNDER, (FP) ; Translate except tone to 
dhs 13e $ mrs _FLOUNDMAT. 
52 53FD 0019 196 TSTG R2 ; test base 
2F 613 QO1C 8 197 BEGL UNDEF JNED 3; undefined 0**0 or 0**(-n) 
54 cf OO1E 198 MNEGL R4, R4 3; R4 = exponent: 
29 13 0021 £199 EQL POWGIX : if exponent is 0, return RO = 1.0 
0093 301 :+ 
0093 20¢ ; Exponent is > 0 or (exponent is =< 0 and base is not = 0 -- use ‘exponent:) 
0 04 ° 
00 £4 8 3 05 phe BBSC #0, R4 peat im, ; branch if ‘exponent: is odd 
8F 9C 0027 06 SQUAR: ROTL al, af 3; R4 =_lexponent:/2 
52 44FD 002C 07 SQUARI : MULG2 ; R2/R3 = current power of bate 
0 8 88 3: Floating overflow will c feuit 
0 : and signal SS$_FLTOVF or 8 sr TOVF LF. 
54 =€9 Ops 1 BLBC R4, SQUAR 3; branch if next bit in Ee esietl is 0 
003 1a 3¢ 
0 13 : Here when bit i of ‘exponent: is a 1. 
13 ; Partial result = partial result * (base * 2**i) 
0 16 ° 
3 19 PARTIAL: 
52 44FD H MULG2 R2, RO ; RO/R1 = new partial result 
BF 18 0 1 A #-1, R4, R4 : R4 = texponent!/2 
FE 1 6 ¢ 30 BNEG SQUAR1 - : loopback if more exponent bits are 1 
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OTSSPOUGJ = G AEATSB se GRVEGERSE grSEpa1o8e 1saBis LATAN TS eesereoboey mans1 P28e (8, 
8 gue helene oettan n 
+4 exponen » return 
6D 6C'AF 9E MOVAB B“EXC. ate exceptions te 


HNDLR_OVER, (FP) ; Trans 
; MTHS_FLOOVEMAT. 
RO, #1, RO : RO/RT = 1.0/result 


DIVG3 
POWGJX: RET ; return, result in RO 


3% 
3; Undefined exponentation error - 0**0 or 0*##(-n) 


UNDEF INED: 


#15, #1, RO 
MOVZBL #MTHSK UNDEXP, =(SP) 
CALLS #1, G*ATHSSSIGNAL 


5 ake lel lated de eee et ak 
VCUUIOTT9FNNNWwem 


50. 01 =6OF 6079 
OO'S8F 9A 


7 RO/R1 = reserved floating operand 
00000000'GF 01 FB 


; Indicate undefined exponentiation. 
3; convert to 32-bit condition code 
: and SIGNAL MTH$_UNDEXP 
3; Note: 2nd arg not needed since no JSB OTS$ 
; is possible. 
RET ; return 
g* 
; The following handler is established to process exceptions which imply 
; underflow of the final result (floating overflow in either of the two MULG's 
3 wen or KA” On the occurrence of such an exception, the handler signals 


EXC_HNDLR_UNDER: 
.Q@ORD “M<R2, R3, R4> : Entry mask 
SETUP ; Set up RO:R3 and identify condition. | 
; Return only if FLTOVF or FLTDIV. 
BBC #PSLSV_FU, SFSW_SAVE_PSW(R2), CON_U 
; Branch Tf caller has not enabled FU. 
poze. +404 1 derma R4 ; Report MTH$_FLOUNDMAT, not SS$_FLTOVF. 


28 
1—E 04 A2 06 €1 
54 00°8F 9A 

oc 


3+ 
3; The peed handler is established to process exceptions which nay & 
; overflow of the final result (floating overflow in either of the two MULG's 
; when exp > 0, floating overflow in the DIVG, or ag bene divide by zero in the 
; DIVG). On the occurrence of such an exception, the handler signals 

3 MTHS_FLOOVEMAT. 


EXC_HNDLR_OVER: 
.QORD Entry mask vce 

Set up RO:R3 and identify condition. 

Return only if FLTOVF or FLIDIV. 

Make the default result -0.0. 

Report MTH$_FLOOVEMAT, not SSS$_FLTxxx. 


; Report caller's PC, not exception PC. 
; Report MTH$_xxx, not SS$_xxx. 

: Signal the Condition. 

; if cont inyed. restore RO and R1. | 

; Continue from the original exception. 
; Exit from handler. 


“M<R2, R3, R4> 
SETUP 


ASHL #15, #1, RO 
MOVZBL #MTHSK_FLOOVEMAT, R4 
DO_SIG: PUSHL — SFSL_SAVE_PC(R2) 
#2, G°MTHSSSIGNAL 


CHF$L_MCH SAVRO(R3) 
S*#SS$_CONTINOE, RO 


SNS AAPA AA AAD TTT TU BB BB BEE BIN NII NI IIInppnonony 
NOUS WN 3S OO NAME WIN OOO NA UE WN 9 ODNAU EWN 0 OO NOU EWR OOWWNOULS WP 
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1-005 OTSSPOWGJ = G REAL*8 ** JNTEGER*4 66-SEP-1984 11:28:16 CMTHRTL.SRCJOTSPOWGJ.MAR;1 
| 7 
j s¢ 
¢ 5 ; Common setup routine for handlers. Returns normally if exception was FLTOVF, 
C ¢ ; FLTOVF_F, FLTDIV, or FLTDIV_F. If the exception was engthing else 
O8C : euequtTs a RET seusio an ex't from the handler with RO = 0, which is 
08C 4 : equivalent to §$ -RESIGNAL. in the case of a normal return (FLTOVF, FLTOVF_F, 
age $3 wLTOLY, or FULTRIVAE? it sets up RO:R3 as follows: 
O8C $ 3 : address of establisher's frame 
008C ss 3 R3: address of mechanism array 
| 09 c $59 
i 
50 7C¢ 5 C 91 SETUP: CLRQ RO ; Set default result to 0.0. 
52 04 AC 7D 008 38 MOva CHFSL_SIGARGLST(AP), R2 ; %2 = address of signal array 
| 009 9 3; 45_= address of mechanism array 
0000'8F 04 A2 B81 009 94 CMPW = CHF SL_SIG_NAME(R2), #SS$_FLTOVF 
04 95 3; Was it a floating overflow trap? 
18 13 0098 96 BEQL DO_RSB ; Branch if yes. 
| 0000'8F 04 A2 Bi 009A 297 CMPW © CHFSL_SIG_NAME(R2), #SSS$_FLTOVF_F 
00A0 98 ; Or a Floating overflow fault? 
10 13 QOA0 99 BEQL DO_RSB ; Branch if yes. 
0000'8F 04 A2 B81 OQOA2 00 CMPW CHF$L_SIG_NAME(R2), #SS$_FLTDIV 
BpAe 01 ; Or a floating divide by zero trap? 
8 13 QOA8 6 BEQL DO_RSB : Branch if yes. 
0000'8F 04 A2 Bi OOAA 30 CMPW © CHFSL_SIG_NAME(R2), #SS$_FLTDIV_F 
ie 04 ; Or a Floating divide by zero fault? 
D9 = «12 Boos ti BNEQ DO_RET ; mere of the above: return from handler 
; Wl = . 
08 A2 97 AF OE 0089 07 DO_RSB: MOVAB B*POWGJX, CHFSL_SIG_NAME+4(R2) 
008 08 ; Change return PC to POWGJX. 
52. 04 A3 DO 00B7 09 MOVL CHFSL_MCH_FRAME(R3), R2 ; R2 = address of establisher's frame 
05 00BB 310 RSB ; Return. 
0oBC 8 311 
ooBC 8312 END 
j 
i 
i 
' 
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“OTSSPOWGJ - G REAL*8 ** INTEGER®4 power routine 1o-SEP-1984 01:59:33 YAX/WMS Macro v04-00 Page 8 ‘2 
| Symbol table 6-SEP=-1984 11:28:16 (CMTHRTL.SRCJOTSPOWGJ.MAR; 1 (4) | 
BAS = 00000004 
CHFSL_MCH_FRAME= 00 4 
CHFSL_MCH_ SAVRO= : ¢ 
CHFSL-SIGARGLST= 000 4 
CHFSL_SIG_NAME = 000 4 
CON_U i464 4 R 0 
DO_RET 0000088 R 0 
DO_RSB 0000 R 0 
DO0_SIG 00000 2 8 
EXC_HNDLR_OVER ae a'eale t8 fe 
Exe HNDLR “UNDER Bepnenae K 0 
= vv | 
EXPGTR 90000025 R 02 
MTHS$$S1GNAL eereeere§ =X 8600 
MTIHSK_FLOOVEMAT *eeeeeee x §° 
MTHSK_FLOUNDMAT eeeeeee x 0 | 
MTHSK UNDE XP eeeerere§ =X 00 
oTs J 00800000 RG 02 
PARTIAL 00000033 R 02 
GJX 00000046C R 02 
PSLSV_FU = 00000006 | 
TuP 0000008C R 02 
SFSL_SAVE_PC = SOO ORAL 
SF SW" SAVE “PSw = 0000 
SQUAR 090000¢7 R 68 
SQUARI 0000002C R 0 
$S$_CONT INUE eterrere x 00 
SSS_FLTIDIV eeeertene x 00 
SSS_FLIDIV_F eeeereee = X 00 
5S$_FLTOVF teeeeeee xX 00 
S$-FLTOVE _F eeeeeeee =X 00 
UNDEF INED 0000004D R 02 
gm er ewe ene mene nm + 
! ;_Psect synopsis ! | 
PSECT name Allocation PSECT No. Attributes | 
» Ob .4 00000000 < 0.) 00 ¢ O.) NOPIC  USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE | 
SABSS$ 00000000 < 0.) O1¢ 1.) NOPIC USR CON’ ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_OTSSCODE ooooocdodBcC ¢ 188.) 02 ¢ 2.) PIC USR CON’ REL LCL SHR EXE RD NOWRT NOVEC LONG 
baw one we mee mee were mm newman $ 
! Performance indicators ! | 
Phase Page faults CPU Time Elapsed Time | 
Initialization 5 00:00:00.09 00:00:00.97 
Command processing 128 $0; 90; 93 00:00:08 3 
Pass 1 1 00:00:02.04 0:00:08.1 
Symbol table sort g 0:00:00.09 0:00:00.16 
Pass 2 6 :00: 36 0:00:03.90 
Syabot cable euteut : 4 :00: ; . : . 8 
sect synopsis outpu :00:00. :00: 
boonmunbterense outoatt 0:00. § 0:00:00 06 


| 
| 
| 
| 
| 
| 
: 
| 
| 
| 
| 
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| VAK=-11 Macro Run Statistics . eee Saggcebgee mie) 7 t 38; i 


Assembler run totals 376 00:00:03.94 00:00:17.21 


The working set Limit was 1050 pages. 
1 bytes (18 pages) of virtual memory were used » buffer the intermediate co 

There were 10 pages of symbol ig space al ocarg? to hold 106 non-local and 8° local symbols. 
source Lines were read in Pass 1, producin $ obs ject records in Pass 2. 

11 pages of virtual memory were ~a to Get ines 10 macros. 


werrreree- erm ee eee eee ewone oe ee 


H ; Macro Library statistics ! 


eee 


Macro library name Macros defined 
“$255$DUA28: CSYSLIBJSTARLET.MLB;2 6 
148 GETS were required to define 6 macros. 


There were no errors, warnings or information messages. 


LSS sss sss sss sss ss nse 


| 
| 
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